New Pseudo-Phase Structure for a-Pu 
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In this paper we propose a new pseudo-phase crystal structure, based on an orthorhombic dis- 
tortion of the diamond structure, for the ground-state a-phase of plutonium. Electronic-structure 
calculations in the generalized-gradient approximation give approximately the same total energy for 
the two structures. Interestingly, our new pseudo-phase structure is the same as the Pu 7-phase 
structure except with very different b/a and c/a ratios. We show how the contraction relative to 
the 7 phase, principally in the z direction, leads to an a-like structure in the [0,1,1] plane. This 
is an important link between two complex structures of plutonium and opens new possibilities for 
exploring the very rich phase diagram of Pu through theoretical calculations. 



Pure plutonium (which also includes Pu with small 
amounts of other solute atoms) is an exotic metal that 
is poorly understood. It appears to be going through a 
Mott-insulator transition in its 5/ electronic orbitals as it 
transforms into different crystal structures as a function 
of temperature. Actinide elements lighter than Pu have 
itinerant 5/ electrons; elements heavier than Pu form a 
second rare-earth series with localized 5/ electrons. The 
a-Pu phase (zero-temperature ground state) is believed 
to have mainly itinerant 5/ electrons. However, 5-phase 
(fee) Pu is believed to be in neither limit (neither com- 
pletely itinerant nor localized), but is instead in some 
kind of exotic electronic state that we may describe as 
partially localized, for want of any better terminology to 
apply to this system 

The understanding of this phase transformation, with 
its associated large volume per atom expansion, has been 
an exciting topic of discussion within band-structure the- 
ory and related correlated-electron extensions to it for 
several decades now. But almost all of the more sophis- 
ticated recent calculations have been done for the fee S 
phased S II 13 

This crystal structure has one atom 
per unit cell and is amenable to better, more computa- 
tionally intensive techniques for handling the correlated 
electron physics, which is believed to be the fundamen- 
tal issue controlling the /-shell Mott-insulator transfor- 
mation. Because the physics is so difficult, we believe 
that studying the strongly correlated 5 phase in isolation 
is less useful than also including the more weakly cor- 
related a phase in the study. A correct theory should 
get both limits right, and is a strong constraint on the 
theoretical approaches that are used. The chief problem 
with including the a structure is its complicated crystal 
structure, which is a low symmetry monoclinic structure 
with space group P21/m and involves 16 atoms per unit 
cell0,0|. Because of this complexity, for example, a fully 
relaxed structure (volume and internal parameters) has 
not yet been obtained with DFT. 



We believe that we have now found a solution to this 
problem. In this paper we propose a new pseudo struc- 
ture for a Pu that involves only two atoms per unit 
cell and is a relatively symmetrical orthorhombic crys- 
tal structure. We show that the real a crystal structure 
only involves small distortions of various atoms around 
this pseudo structure and that the electronic-structure 
energy of both systems (at least at the band-structure 
level) is nearly identical. 

By a pseudo structure we mean a simpler crystal struc- 
ture where the atoms are quite close to the positions of 
the more complex real structure. Pseudo-structures are 
interesting for two reasons: (1) they provide insight into 
the geometrical structure of the atoms in the original 
structure and possible phase-transformation pathways, 
and (2) they provide an alternative less expensive crys- 
tal structure for complicated electronic and microstruc- 
tural calculations. From a physics point of view, one can 
view the pseudo structure as an approximation to the 
real structure that captures the dominant total energy 
of the electronic-structure (heat of formation) and local 
bonding effects; the real structure can then be viewed as 
just a minor distortion of the pseudo structure, with only 
a small energy difference between the two structures. 

In the case of Pu, insight into the crystal structure is 
extremely important. Unlike, for example, the cubic di- 
amond structure of Si where it is easy to understand the 
tetrahedral nearest-neighbor structure as a consequence 
of directional s-p bonds, or the close packed fee (face- 
centered-cubic) crystal structures favored by Madelung 
energies in metals, it is extremely difficult to under- 
stand why the lowest temperature Pu phases are so ex- 
tremely anisotropic and complex. In a Pu, for example, 
the nearest-neighbor structure involves a complicated ar- 
rangement of short and long bondsQ- How can we un- 
derstand this, and why do the crystal structures become 
more isotropic and metallic like at high temperatures? 
Good pseudo structures give us the opportunity to look 
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at these crystal structures in new ways that may eventu- 
ally stimulate insights that will resolve some of the ques- 
tions. 

Understanding the complex phase diagram of Pu is 
likely to involve consideration of strong directional bond- 
ing and high density of states effects of the / electrons 
and how that bonding weakens under strong electron- 
electron correlations. Sophisticated electronic-structure 
calculations are necessary to make progress in this area. 
Unfortunately, the a and phases, with 16 and 34 atoms 
in the unit cell, respectively, are just too computationally 
expensive for such calculations to be practical. To have 
a simple pseudo-phase crystal structure that can mimic 
the actual bonding arrangement of the extremely compli- 
cated crystal structure of the real phase makes it possible 
to do these very sophisticated electronic-structure calcu- 
lations with little loss of accuracy or insight. If accurate 
pseudo phases exist, they also make possible complicated 
microstructural models of the phase transformations that 
captures the dominant components of free energy differ- 
ences without requiring huge numbers of parameters. For 
example, Ginzburg-Landau models become much more 
accessible through the use of pseudo phases. An attempt 
is currently being made to do this for the very important 
a to 5 phase transformation of Pu0- 

What is even more interesting is that our proposed 
pseudo structure is identical with the 7 phase of Pu, ex- 
cept that the unit cell dimensions are greatly changed. 
This crystal structure is an orthorhombic distortion of 
the cubic diamond lattice. It has two atoms in the small- 
est unit cell, one at the origin, and the second 1/4 of the 
way along the a, b, and c axes (which are now all differ- 
ent in the orthorhomic unit cell). Because both atoms 
have a face-centered configuration, there are 8 atoms in 
the primitive orthorhombic unit cell (just like cubic dia- 
mond). As we will show, our pseudo structure is far bet- 
ter than the only other previously proposed pseudo struc- 
ture for a Pu, the Crocker pseudo structure^. This 
structure was based on the insight that the the a struc- 
ture can be viewed as two repeating planes of a distorted 
hexagonal structure. Although this structure also has 
two atoms in its unit cell (each of the distorted hexagons 
are replaced by symmetric hexagons), it is also mono- 
clinic, and has not proved useful for electronic structure 
calculations. We will also show that its near-neighbor 
shells are a poor approximation to the real a phase. 

What makes Pu fascinating from an electronic- 
structure point of view is the instability in the 5/ shell, 
which have itinerant electrons for a Pu, and localize into 
the core for Am and higher-Z actinides. This Mott- 
like transition for the 5/ electrons involves very strong 
electron-electron correlations, and requires a much more 
complicated theory for the electrons than the usual lo- 
cal density approximation (LDA) for density functional 
theories (DFT), which are like a mean-field treatment of 
the correlations llj. If the parabolic decrease of the vol- 



ume of the light actinides as a function of atomic number 
can be explained with a simple Friedel model, plutonium 
deviates from this theoretical estimation, since Pu has 
a larger volume than Np. Until now there is no defini- 
tive explanation of this upturn. In addition, it is well 
known that DFT in its local density or generalized gradi- 
ent approximation (GGA), due to the presence of strong 
correlation effects, has been unable to describe the 5- 
phase of plutonium, underestimating its atomic volume 
by about 20%, for example. But recently, several at- 
temps to go beyond LDA, have given some understand- 
ing of the a-S transition. The LDA+C/ method0,3 and 
the LDA+DMFT (dynamical mean-field theory) [3j, by 
taking into account the localization of the /-states, have 
shown the role of spin and orbital polarization in describ- 
ing the (5-phasc. 

Soderlind et al\l^ have performed ab initio relativistic 
total energy calculations and offered an explanation for 
the crystal structure of a-Pu in terms of a Peierls-Jahn- 
Teller-like distortion. Comparing the total energy of the 
a phase with total energies of symmetric (fee, bee, hep) 
structures and other actinides- like structures (such as a- 
U and a-Np) they showed that the a-Pu structure min- 
imizes the total energy for plutonium, but they couldn't 
explain the upturn in equilibrium volume for Pu. More 
recently Jones et q^|l3| have performed calculations with 
two different methods, FLAPW and LCGTO-FF for all 
the light actinides. Although their theoretical volumes 
agree well (within about 2%) with the experimental vol- 
umes for Th to Np, they show a larger discrepancy (7%) 
for a-Pu, and no upturn between Np and Pu. No at- 
tempt was made to relax the structural parameters of 
the a phase. They also speculated that some localiza- 
tion effects could already be present in a-Pu accounting 
for the discrepancy between the computed and measured 
atomic volume and causing the upturn. To our knowl- 
edge, no fully relaxed electronic structure calculations 
have been performed on the a, (3 and the 7 phase and 
a complete and comprehensible picture of the plutonium 
phase diagram still does not exist. 

In our calculations, we have mainly used the projector 
augmented-wave (PAW) method, which combines ideas 
from pseudopotentials and the linearized augmented- 
plane-wave (LAPW) method ^ij. Forces and the full 
stress tensor can be easily calculated and used to re- 
lax atoms into their instantaneous ground state. This 
is implemented in the plane-wave basis VASP code [3 
with spin-polarized DFT and the GGA of Perdew and 
Wang|16j. Ultrasoft Vanderbilt pseudopotentials were 
used to represent the core electrons. 

This allowed us to perform for the first time a fully 
relaxed calculation of the a-phase. Results are presented 
111 Tabled As in previous calculations E3- the atomic 
volume is slightly lower compared to the experimental 
value. This discrepancy is partially due to the absence 
of spin-orbit coupling in our calculations 0, 0|. The 
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relaxed structure is very closed to the experimental one 
(the energy difference is about 4 mRy), and the positions 
of the atoms relaxed by only 2 or 3 %, showing that DFT- 
GGA reproduces the ground state of plutonium with a 
good precision. 

We have compared the total energies of the a structure, 
Crocker's structure, and the 7 structure. As expected, 
the a phase has the lowest energy, while Crocker's struc- 
ture is about 20 mRy larger, and the 7 phase 36 mRy 
larger. The theoretical volume of the 7 phase is very 
far from the experimental volume (23.14 A 3 ), but very 
close to the experimental volume of the a phase, see Ta- 
ble n If we minimize the total energy of this structure as 
a function of b/a and c/a we obtain a large contraction 
in the z direction and a smaller one in the y direction 
to obtain a structure very different from 7 (we will call 
this new structure "pseudo-a" in the rest of the paper). 
The volume increases slightly and is roughly equal to 
the theoretical volume for the a structure, see Table [i] 
Surprisingly the total energy is also very close to the to- 
tal energy obtained for the a structure, about 3 mRy 
higher, almost within the accuracy of calculations. Us- 
ing a full-potential linear augmented-plane wave method 
(FLAPW) we have checked the influence of the spin-orbit 
coupling on our results (see TableHJ). The total energy of 
the pseudo-structure is still very close to the a one and 
the spin-orbit coupling doesn't help to differentiate the 
two structures. 

It is interesting to speculate on why we find nearly the 
same total energy for two structures that appear to be 
so completely different. To understand this it is useful to 
compare the interatomic distances in the 7, the pseudo- 
a, and the a structure, which are shown in Table [H] 
Although the atomic environment around each atom in 
the a structure is very complex, we can divided it into 
roughly two categories: short and long bonds, well sepa- 
rated. Each atom has four first-nearest neighbors and ten 
second-nearest neighbors, except atom I which has 3 and 
13 and atom VIII which has 5 and 7. During the trans- 
formation from 7 to pseudo-a the first shell of neighbors 
contracts around the atom while the second and the third 
shell expand to be very close to the fourth shell. This 
gives a local environment of 4 nearest neighbors and 10 
second-nearest neighbors for the pseudo-a , exactly what 
is found in average in the a structure. In addition, the 
interatomic distances are very close. So the relaxation of 
the 7 phase gives us a new structure, with nearly exactly 
the same distribution of interatomic distances than in the 
a phase. Crocker's structure, however, is quite different 
with only two first-nearest neighbors, and an intermedi- 
ate shell with four atoms. So, the local atomic environ- 
ment is more important than the monoclinic symmetry 
in determining the energy of the a structure. 

In Fig. 01 we show the density of states (DOS) ob- 
tained for the three structures at their theoretical vol- 
umes. They are very similar, with the Fermi level in a 



pseudo gap and two broad features, one above and one 
below the Fermi level. The DOS of the 7-Pu with the ex- 
perimental parameters has been compared to the pseudo- 
a structure and we can see the presence of a peak at the 
Fermi energy, so during the contraction of the cell along 
the b and c axis there is the opening of a large pseudo gap 
with the reduction of the density of states at the Fermi 
level. 

Looking at the [Oil] plane of the pseudo-a structure, 
see Fig. |2J one can almost visualize the a phase. The 
atoms in dark and grey are the two distinct atoms in the 
pseudo-a structure, and we have superimposed in blue 
the a-structure atoms. They are very similar despite 
the complexity of the a pattern and the relative simplic- 
ity of the distorted diamond structure. We have calcu- 
lated the distance between related atoms in the two dif- 
ferent crystal structures, which has an average of about 
0.25 A. These distortions are responsible for the small 
total-energy difference between the two structures. In 
addition, the [Oil] planes have the same interlayer dis- 
tance as in the a phase (2.31 A and 2.26 A for pseudo-a 
and a respectively). As we indicated earlier the new 
structure respects the average surrounding of an atom in 
the a phase, with four nearest neighbours, but exactly 
as in the a phase, there is two in the same plane, one in 
the plane above and one in the plane below. The angle 
between one atom and its nearest neighbours from the 
two other planes is also very similar (around 137° and 
130° for pseudo-a and a respectively, in diamond it is 
109°). In addition, the connection of our pseudo phase 
with the gamma phase suggests a deep structural rela- 
tionship between the low-temperature phases of Pu that 
needs further examination. It is interesting to observe 
that the b/a and c/a ratios of the real 7 phase of Pu 
have a tendency to move toward the values for our new 
pseudo phase as the temperature is lowered. Zachariasen 
and Ellinger measured the lattice constants in function of 
temperature. The lengths of b and c both decrease with 
temperature, while the length of a increases^. If the (3 
phase did not interrupt this trend, it is likely that 7 Pu 
would slowly continue in this way until it would reach the 
phase transformation into the a phase that would then 
require just a tiny distortion. 

In conclusion we have discovered a very simple pseudo 
structure reproducing the main features of the complex 
ground state of plutonium. Our results show that the 
local atomic environment in the a structure with four 
first-nearest neighbors and 10 second-nearest neighbors 
is its crucial characteristic. We have found a very simple 
path (following the experimental variation of the lattice 
constants of 7 Pu with temperature) between 7 and a 
plutonium. These results show the predictive power of 
DFT calculations and our new pseudo structure of the a 
phase opens new possibilities for exploring the behavior 
of plutonium. It will allow approximate calculations of 
the a phase with far more sophisticated techniques than 



4 



TABLE I: Structural and internal parameters (A) for the ex- 
perimental and theoretical a, Crocker, pseudo-a and 7 struc- 
tures of plutonium. For a, the volume and all the internal 
parameters have been relaxed. The energies (mRy) are rela- 
tive to the a structure, which is the ground state for all the 
calculations, with or without spin-orbit coupling. 



Pu-a Crocker 7 pseudo-a 





Exp. 


Th. 








a 


6.19 


5.96 


2.94 


2.85 


3.50 


b/a 


0.78 


0.76 


1.56 


1.83 


1.61 


c/a 


1.77 


1.83 


1.93 


3.22 


2.09 




101.8 


101.8 


116.3 






Volume 


20.0 


18.0 


17.3 


17.0 


18.0 


Energy (CCA) 







20 


36 


3 


Energy (GGA+SO) 







18 
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TABLE II: Nearest neighbours distances (A) for the Crocker 
structure, the 7-Pu, the pseudo-a, and the a-Pu. These val- 
ues are given for the theoretical equilibrium volume for each 
structure. 

distances (number) 
Crocker 7-Pu pseudo a 

1 n.n 2.44(2) 2.73(4) 2.47(4) 2.46-2.55(4) short bonds 

2 n.n 2.79(2) 2.85(2) 3.31(4) 

3 n.n 2.79(2) 2.97(4) 3.50(2) 3.21-3.59(10) long bonds 

4 n.n 3.02(2) 3.39(4) 3.49(4) 



have been previously possible. 
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FIG. 1: Density of states for a-Pu (upper pannel), pseudo-a 
and 7-Pu in dashed (middle panel), and Crocker's structure 
(lower panel) calculated at the theroretical equilibrium vol- 
ume. 
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FIG. 2: Distorted diamond structure (pseudo-a) in the (011) 
direction. The blue atoms are the eight positions in the a 
structure. 
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